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, Abstract 

i~{ . It is shown how the phase-damping master equation, either in Marko- 

• vian and nonMarkovian regimes, can be obtained as an averaged random 

unitary evolution. This, apart from offering a common mathematical 
00 , setup for both regimes, enables us to solve this equation in a straight- 

forward manner just by solving the Schrodinger equation and taking the 
stochastic expectation value of its solutions after an adequate modifica- 
tion. Using the linear entropy as a figure of merit (basically the loss 
of quantum coherence) the distinction of four kinds of environments is 
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. suggested. 

O ■ 1 Introduction 

I The principle of quantum superposition is at the base both of important issues 

■ concerning the conceptual foundations of quantum mechanics and of poten- 

CIhI tial direct technological applications of quantum theory |§, ^ . In particular the 

preservation and dissapearance of these quantum superpositions play a central 
C I role in topics like the quantum-to-classical transition Q and the performance 

of some quantum information processing tasks |^ . Physically this decoherence 
has its origins in the pervasive effect of the surrounding environment. Mathe- 
matically it is expressed by means of master equations (ME hereafter), that is 
evolution equations for the reduced density operator of the system which are 
usually obtained after some careful assumptions M. 



By and large two main approaches can be followed to study these master 
equations namely a constructive and an axiomatic approach. In the former 
one begins considering the joint system-environment degrees of freedom and 
finally pursues reasonable and manageable approximations to obtain the master 
equation of the system. This approach usually falls into mathematical difficulties 
hard to deal with, though some limit situations have been satisfactorily dealt 
with 1^, |j. On the other hand, the latter approach poses a set of reasonable 
and physically motivated assumptions |l^ and then derive the general form 
for the master equation of systems satisfying that set of axioms. Though most 
of ME of practical relevance fall into this category (see e.g. [Q), the range of 
applicability of this approach is limited by definition, especially regarding one 
of its fundamental axioms, namely the semigroup condition (Markovianity) . 

In an attempt to embrace also nonMarkovian situations we have combined 
well-known stochastic methods jointly with common tools of functional analy- 



sis employed in quantum mechanics |13, to state the following result: Any 
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Lindblad-type evolution with selfadjoint Lindhlad operators, whether Marko- 
vian or nonMarkovian, can be understood as an averaged random unitary 
evolution. This result, whose fully-fledged mathematical proof is provided else- 
where ||l5|, reveals the fact that the Lindblad structure of ME for the case of 
selfadjoint Lindblad operators is beyond the Markov approximation and can 
also be recovered in nonMarkovian situations. 



Here we exploit the mathematical advantages of this approach to study the 
loss of coherence described by the phase-damping ME and identify different 
types of environments using the linear entropy as a figure of merit. 

2 The Phase-Damping Master Equation as an 
Averaged Random Evolution 

The generalized phase-damping ME reads 



^ = -^[H,p{t)]-^[H,[H,pm (1) 

where H denotes the Hamiltonian of the system and we have allowed X{t) to 
be any nonnegative real function. In general eq. (Q) denotes a nonMarkovian 
evolution reducing to a Markovian one only when A(<) = A = const. Eq. (||) 
is a Lindblad-type ME with selfadjoint Lindblad operators and consequently 
following the result stated above can be obtained as an averaged random unitary 
evolution. The proof for this case is straighforward. Let U{t) = e^**-^ (fi, = 1) be 
the unitary evolution operator of the quantum system under study. Let us now 
define the random unitary evolution operator Ust(t) — e~^*^~^^ ^ where 

a{t) denotes an arbitrary real function and Bt denotes standard real Brownian 
motion ||l6|, thus a{s)dBs will denote the well-known Ito integration. Now 
let the density operator p{t) be defined as 



p{t)^E[Ust{t)pmUt)] (2) 

where E denotes expectation value with respect to the probability measure of 
Bt- Then after making some calculations one can see that p{t) satisfies the ME 
(0) with X{t) = a^{t). 



3 Interference Damping for Plane Waves 

One of the immediate consequences of the above procedure arises by realizing 
that the solutions to the ME (|l|) can be obtained by performing the substitution 
i — !■ t + a{s)dBs in the unitary solutions and then taking the stochastic 
expectation value. 

We will illustrate these ideas with the study of the interference pattern be- 
tween two equally- weighted plane waves with momentum fci and ^2 respectively 
evolving under the ME (M. The unitary solution can be obtained elementarily 
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either by solving the Schrodinger equation for the state vector or by solving the 
Liouville-von Neumann equation for the density operator or by resorting to the 
Wigner characteristic function xi-^j — Tr[/9(t)e*'-'^P"'"'''''^] where q and p are 
the position and momentum operators (cf. also The interference pattern 

is provided by the well-known expression 



^2 _ ^2 

lix, t) = 1 + cos[(fci - k2)x i- -t] (3) 

2m 

This pattern is changed when the initial superposition evolves under eq. 
and accordingly the solution is more involved. One can again resort to the 
Wigner characteristic function x(A, /i) to obtain the solution iQ, but using our 
previous result one can immediately arrive at the desired expression: 

I{x,t) ^ i{x,t) =E[I{x,t+ f a{s)dBs)] (4) 

Jo 

Then knowing that (cf. appendix) 
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(5b) 



we get 



-1 f^l^V f>'Mds k?-k1 

I{x,t) = l + e '"^ J cos[(fci ^ fc2)a; - \^ A (6) 

Different interference patterns are obtained corresponding to different types 
of environments. Note that since the Markovian approximation is not assumed, 
we have now gained some generality. 



4 Different Types of Environment 

We will investigate the time evolution of the linear entropy in order to study 
the progressive loss of coherence of systems driven by the ME (0). To avoid 
mathematical nuisances because of the lack of normalizability of plane waves, 
let us focus upon an initial gaussian wave packet 

ij{x,0) = 



(2^a2)i/4 

Under equation (0) the linear entropy Sun (t) of the system can be calculated: 
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1 4TOcr^ dy 



= -fe (7) 

where X{t) = J* a'^{s)ds. Note that the previous integral exists for all possible 
choices of a{t) since the integrand (= g{y,t)) is continuous in M+ and 
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linig(y,<) = —= 
lim g{y,t) = 



giy,t) ~ 2;>1 



Also notice that the second term in (0) (= I{t)) is bounded from above and 
from below: 



0</(t)<1^0<^H„W<l (8) 

In particular Siin{t) = if and only if \{t) = 0, i.e. one always obtains 
a decohering process. Immediately one can distinguish four different kinds of 
situations: 

a) \{t) — At. This is the usual Markovian approximation. In this case (cf. 
also I) 



Sunit) *^ 1 (9) 

b) X{t) Aq. Let us call it SubMarkovian regime. Now the linear entropy 
reaches a maximum value 5*0 < 1, thus the initial coherence is partially 
preserved. 

c) A(t) oo but ^zl^ 0. Let us call it type-I SuperMarkovian regime. 
Now as in the Markovian case Siin{t) 1 but at a slower pace. This allows 
the system to preserve a higher degree of coherence at a given time. 

d) X{t) * — °? oo and t^l^^ qo. Let us call it type-II SuperMarkovian 
regime. This is the most destructive case since Sun{t) 1, that is co- 
herence is completely lost and even faster than in the usual Markovian 
case. 

The four situations are qualitatively represented in Fig. 0. 
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Figure 1: Decoherence functions for different types of environment. See text. 



5 Conclusions 

We have derived the phase-damping master equation using stochastic calcu- 
lus and subsequently solve it using the advantages of this procedure. These 
advantages run from embracing under the same mathematical formalism both 
Markovian and nonMarkovian evolutions (thus allowing us to deal with more 
general physical situations) to providing a new method to solve this equation. 
The great advantage of this procedure comes out from the fact that you only 
need to solve the unitary evolution equation, which can be done for the state 
vector instead of dealing with the density operator. This represents a compu- 
tational advantage. Once the unitary solution is known all one has to do is to 
find the stochastic expectation value of this solution after modifying it with the 
appropiate random factor (see above). 

Using then the linear entropy as a figure of merit we have found four types 
of environments according to the rate of loss of coherence they induce upon 
the system. In particular it may be concluded that what we have called Sub- 
Markovian and Type-I SuperMarkovian regimes appear as better choices to per- 
form quantum-informational and quantum-computational tasks than the usual 
Markovian environment. Now the physical conditions should be provided to 
properly identify what these kinds of environment are. 

It should be noticed that as a by-product of this approach, independently of 
the function (T(t) used to describe the system-environment response, a positive 
and nondecreasing rate of decoherence is always obtained, which indicates that 
it is an irreversible process. Thus under the conditions in which the phase- 
damping master equation is valid to describe the evolution of a system, coher- 
ence can never be recovered from the environment, whether this is Markovian 
or nonMarkovian. 
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A Some Computational Notes 



In order to compute the expectation values ( pa| ) and (5b) it is necesarry to know 

the nth order moments of the real- valued stochastic process 

It is simple to find /?„(<) = E[Xj"] using Ito's formula 111 upon (j){Xt) = X^: 



- n/3„_i(t) + ^^i^^/5„_2(iV'(i) n>2 (10) 
Since /3o (t) = 1 and /3i (t) = it is straighforward to recursively arrive at 



p2n{t) - |^A"(t) (11a) 

f32n+l{t) = (lib) 

where A (t) = J^a^{s)ds. 

To arrive at (|^) the following integral is needed which can be obtained by 
means of already tabulated integrals [M : 



e 



d0=J-CTi{^) (12) 
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